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Abstract 


Progress  is  reported  in  development-,  implementation,  and  application  of  a  spectral  method 
for  ab  initio  studies  of  the  electronic  structure  of  matter.  In  this  approach,  antisymme¬ 
try  restrictions  are  enforced  subsequent  to  construction  of  the  many-electron  Hamiltonian 
matrix  for  an  atom  or  molecule  in  an  orthonormal  spectral-product  basis.  Transformation 
to  a  permutation-symmetry  representation  obtained  from  the  eigenstates  of  the  aggregate 
electron  antisymmetrizer  is  seen  to  enforce  the  requirements  of  the  Pauli  principle,  and  to 
eliminate  the  unphysical  (non-Pauli)  states  spanned  by  the  product  representation.  Re¬ 
sults  identical  with  conventional  use  of  prior  antisymmetrization  of  configurational  state 
functions  in  variational  calculations  are  obtained  in  applications  to  many-electron  atoms, 
providing  some  degree  of  confidence  in  the  soundness  of  the  method  more  generally.  In 
applications  to  polyatomic  molecules,  the  development  is  seen  to  provide  certain  potential 
advantages  over  conventional  methods,  and,  particularly,  to  accommodate  the  incorpora¬ 
tion  of  fragment  information  in  the  form  of  Hermitian  matrix  representatives  of  atomic 
and  diatomic  operators  which  include  the  non-local  effects  of  overall  electron  antisymme¬ 
try.  Exact  atomic-pair  representations  of  polyatomic  Hamiltonian  matrices  are  obtained  in 
this  way  which  avoid  the  ambiguities  of  semi-empirical  fragment-based  methods  previously 
described  for  electronic  structure  calculations.  Illustrative  applications  to  the  well-known 
low-lying  doublet  states  in  the  H3  molecule  demonstrate  that  the  eigensurfaces  of  the  an¬ 
tisymmetrizer  can  anticipate  the  structures  of  the  more  familiar  energy  surfaces,  including 
the  seams  of  intersections  common  in  high-symmetry  molecular  geometries.  The  calculated 
H3  energy  surfaces  are  found  to  be  in  good  agreement  with  corresponding  valence-bond 
results,  and  in  general  accord  with  accurate  values  obtained  employing  conventional  high- 
level  computational-chemistry  procedures.  The  spectral-product  approach  to  molecular 
electronic  structure  avoids  the  repeated  evaluations  of  the  one-  and  two-electron  integrals 
required  in  construction  of  polyatomic  Hamiltonian  matrices  in  the  antisymmetric  basis 
states  commonly  employed  in  conventional  calculations  of  adiabatic  potential  energy  sur¬ 
faces,  providing  an  alternative  ab  initio  formalism  potentially  suitable  for  computational 
applications  more  generally. 
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1.0  Introduction 


Considerable  progress  lias  been  reported  over  the  past  four  decades  in  ab  initio  studies  of 
the  complex  (Born-Oppenheimer)  potential  energy  surfaces  which  describe  the  ground  and 
low-lying  excited  electronic  states  of  molecules  [1].  The  contributions  of  Professor  Mark 
S.  Gordon  to  this  enterprise  have  been  remarkably  comprehensive,  as  typified  by  his  role 
in  the  development,  wide-spread  distribution,  and  continuing  collaborative  refinement  of 
the  well-known  GAMESS  suite  of  computer  codes  [2],  Accurate  ab  initio  methods  pro¬ 
vided  in  this  and  other  computational  resources  for  potential  energy  surfaces  commonly 
entail  repeated  calculations  of  large  numbers  of  one-and  two-electron  integrals  over  explic¬ 
itly  antisymmetrized  basis  states  in  construction  of  a  many-electron  Hamiltonian  matrix 
for  a  range  of  atomic  spatial  arrangements,  followed  by  determinations  of  the  energies  and 
associated  eigenfunctions  of  selected  electronic  states  at  each  atomic  arrangement  and  con¬ 
struction  of  the  expectation  values  corresponding  to  physical  properties  of  interest.  New 
methods  for  these  purposes  devised  in  a  spirit  of  continuing  collaborative  improvements 
would  clearly  be  welcome,  particularly  if  they  could  avoid  repeated  calculations  of  the  in¬ 
dividual  electronic  integrals  and  related  quantities  generally  required  in  the  construction  of 
potential  energy  surfaces,  or  possibly  circumvent  entirely  determinations  of  total  molecular 
energies  (and  their  associated  differencing  problems)  in  favor  of  calculations  of  atomic  and 
interaction  energies,  while  still  proving  applicable  to  both  ground  and  electronically-excited 
states  on  a  common  basis. 

In  a  contribution  to  the  systematic  improvement  of  standard  quantum- chemical  methods, 
an  alternative  perspective  is  provided  in  the  present  report  on  enforcement  of  the  antisym¬ 
metry  requirement  on  proper  atomic  and  molecular  wave  functions  [3],  which  requirement 
is  commonly  satisfied  by  the  aforementioned  prior  constraint  on  the  representational  ba¬ 
sis  states  employed  in  variational  calculations.  In  the  present  approach,  rather,  Pauli’s 
exclusion  principle  is  ignored  at  the  outset,  and  antisymmetry  is  enforced  subsequent  to 
construction  of  the  Hamiltonian  matrix  in  a  spectral-product  representational  basis  famil¬ 
iar  from  the  theory  of  long-range  interactions  [4],  When  applied  to  the  electronic  structures 
of  atoms,  the  development  is  seen  to  be  equivalent  to  the  familiar  Slater  approach  adopt- 
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ing  prior  antisymmetry  of  many-electron  configurational  states  functions  [5],  providing 
some  degree  of  confidence  in  the  alternative  method  more  generally.  In  applications  to 
polyatomic  molecules,  the  approach  provides  a  number  of  advantages  over  conventional 
methods,  and,  in  particular,  accommodates  the  incorporation  of  ab  initio  atomic  and  di¬ 
atomic  information  in  Hamiltonian  matrices  in  a  manner  which  avoids  certain  difficulties 
encountered  in  previously  described  semi-empirical  fragment-based  approaches  [6,7]. 

The  present  method  is  ultimately  based  on  conventional  variational  calculations  in  L2  rep¬ 
resentations  of  many-electron  states,  guaranteeing  its  convergence  when  closure  is  achieved, 
in  which  limit  the  distinction  between  use  of  simple-product  and  explicitly  antisymmetrized 
aggregate  basis  states  is  shown  to  become  inconsequential.  As  presented  here,  the  theory 
is  applicable  to  many-electron  atoms  and  to  polyatomic  molecules  which  dissociate  into 
neutral  species  on  their  lowest-lying  potential  energy  surfaces,  but  which  can  otherwise 
involve  arbitrary  admixtures  of  covalent,  ionic,  van  der  Waals,  and  metallic  interatomic 
interactions.  Applications  of  the  approach  to  ionic  systems,  or  to  neutral  systems  which 
can  give  rise  to  ion-pair  states  in  asymptotic  separation  limits,  require  only  minor,  largely 
self-evident,  modifications  of  the  present  description. 

The  atomic  spectral-product  basis  employed  in  the  development  is  known  to  span  the  to¬ 
tally  antisymmetric  irreducible  representation  of  the  aggregate  electron  symmetric  group 
once  and  only  once,  making  it  suitable  for  analytical  and  computational  studies  of  the 
electronic  structures  of  many-electron  systems  [8-12],  The  matrix  representative  of  the 
antisymmetrizer  constructed  in  the  spectral-product  basis  is  employed  in  separating  the 
totally  antisymmetric  and  non-totally-antisymmetric  (non-Pauli)  states  spanned  by  the 
basis,  and  in  correspondingly  isolating  the  physically  significant  block  of  the  Hamilto¬ 
nian  matrix  by  unitary  transformation.  Hermitian  matrix  representations  of  atomic  and 
atomic-pair  operators  are  devised  in  this  way  for  polyatomic  molecules  which  individu¬ 
ally  have  well-defined  asymptotic  separation  limits,  and  which  need  be  constructed  only 
once  for  multiple  applications.  These  attributes  of  the  approach  facilitate  development  of 
an  ab  initio  unitary  transformation  formalism  which  provides  an  exact  atomic-pair  repre¬ 
sentation  of  polyatomic  Hamiltonian  matrices,  including  particularly  modifications  of  the 
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bare  atomic-pair  interactions  due  to  aggregate  electron  antisymmetry  consequent  of  their 
incorporation  in  the  polyatomic  system. 

The  general  theory  is  described  in  Section  2,  where  the  spectral-product  representation 
of  electrons  is  defined,  the  unitary  transformation  formalism  for  isolation  of  the  physical 
portion  of  the  Hamiltonian  spectrum  is  described  for  both  atoms  and  molecules,  and  the 
atomic-pair  form  of  the  polyatomic  Hamiltonian  matrix  is  reported.  Illustrative  calcula¬ 
tions  of  the  familar  low-lying  doublet  electronic  states  of  the  H3  molecule  are  provided  in 
Section  3,  and  concluding  remarks  made  in  Section  4. 

2.0  Theory 

The  spectral-product  approach  to  the  electronic  structures  of  atoms  is  described  in  Section 
2.1  as  an  introduction  to  the  formalism,  its  application  to  polyatomic  molecules  is  presented 
in  Section  2.2,  and  the  atomic-pair-based  implementation  of  the  approach  is  reported  in 
Section  2.3. 


2.1  Spectral-Product  Formalism  for  Atomic  Structure 

Solutions  of  the  Schrodinger  equation  [13] 

H(r)^f(r)  =  ¥(r)  •  E  (1) 


for  an  n-electron  atom  with  Hamiltonian  operator  [14] 


^  i~  2  r-7  2  ^  ^ 

E  i 

i= 1  7=i+] 


(2) 


j=i+ 1  1 ^ 

can  be  constructed  employing  a  formally  complete  square-integrable  (L2)  representation 
in  the  form  [8] 

®(r)  =  {</>(i)  <8>  0(a)  <g>  •  •  ■  0(n)}o  .  (3) 


The  row  vector  comprises  a  denumerable  Hnite  or  infinite  basis  set  of  orthonormal 
one-electron  spin-orbitals  specified  by  the  usual  quantum  numbers  ( n,l,mi,s,ma )  [14], 
where  i  refers  to  the  spin  and  space  coordinates  of  the  ith  electron  and  r  =  (1,  2,...  n) 
represents  collectively  the  coordinates  of  all  n  atomic  electrons.  The  particular  choice  of 
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spin-orbitals  is  irrelvant  for  the  present  development,  so  long  as  the  outer-product  ((g)) 
Hilbert  space  of  Eq.  (3)  can  include  a  suitable  closure  in  the  limit  of  a  complete  spectral 
basis  [15].  The  subscript  “O”  in  Eq.  (3)  indicates  the  adoption  of  a  particular  ordering 
convention  for  the  sequence  of  product  functions  in  the  row  vector  3>(r),  the  consequences 
of  which  convention  are  indicated  when  appropriate  in  the  sequel. 

The  Hamiltonian  matrix  representative  of  the  operator  of  Eq.  (2)  in  the  basis  of  Eq.  (3), 

n  n 

n=  (&(r)\H(r)\&(r))  =  ^  V^j,  (4) 

i= 1  j=i+ 1 

includes  kinetic  =  (#(r)|  —  (h2/2m)V2 |3>(r)),  potential  =  (#(r)|Ze2/rj|3>(r)), 
and  electron-interaction  =  (<F(r)|e2/rjj|<J?(r))  Hermitian  matrix  representatives  of 

the  corresponding  operators.  The  dimensions  of  these  matrices  are  determined  by  that  of 
the  many-electron  basis  of  Eq.  (3),  although  only  standard  one-and  two-electron  integrals, 
and  their  products  with  the  unit  matrices  expressing  the  orthonormality  of  the  remaining 
n-1  and  n- 2  orbitals,  respectively,  are  required  in  evaluations  of  the  Hamiltonian  matrix 
of  Eq.  (4).  The  specific  forms  of  these  matrices  are  determined  by  the  particular  ordering 
convention  adopted  in  Eq.  (3)  [8]. 

The  matrix  Schrodinger  equation  [13] 


H  UH  =  UH  E  (6) 

provides  approximate  or  exact  representations  of  the  row  vector  of  eigenfunctions  of  Eq. 
(1)  in  the  form  ^(r)  =  <&(r)  •  Uh,  as  well  as  the  diagonal  matrix  of  associated  eigenvalues 
E.  The  eigenspectrum  so  obtained  includes  both  physically  significant  (totally  antisym¬ 
metric)  and  unphysical  (non-totally  antisymmetric  or  “non-Pauli”)  solutions,  all  of  which 
are  spanned  by  the  representation  of  Eq.  (3)  [9].  Accordingly,  this  representation,  and 
the  corresponding  Hamiltonian  matrix  and  Schrodinger  equation  of  Eqs.  (4)  and  (5),  can 
encompass  descriptions  of  the  electrons  as  either  Bose-Einstein,  Fermi-Dirac,  or  possibly 
other  forms  of  identical  particles,  necessitating  some  means  of  distinguishing  among  these 
possibilities  [3]. 
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The  physical  and  non-Pauli  eigenstates  obtained  from  Eqs.  (1)  to  (6)  can  be  separated  by 
constructing  a  transformed  Hamiltonian  matrix  which  has  a  physical  block,  a  non-Pauli 
block,  and  a  vanishing  off-diagonal  block.  This  separation  is  accomplished  by  constructing 
the  eigenstates  of  the  n-electron  antisymmetrizer  in  the  representation  of  Eq.  (3)  and 
transforming  the  Hamiltonian  matrix  of  Eq.  (4)  to  this  new  representation.  The  required 
representation  matrix  is 

PA=(*(r)|P.4l*(r)>,  (7) 

where  the  Hermitian  antisymmetrizer  employed  [3], 

n! 

Pa  =  El-T'h,  (8) 

P=  1 

is  left  unnormalized.  The  matrix  of  Eq.  (7)  is  found  to  comprise  a  series  of  individual 
non-zero  n!-by-n!  blocks  on  the  diagonal  the  elements  of  which  are  all  either  +1  or  -1,  and 
a  remaining  set  of  vanishing  blocks.  The  dimensions  and  natures  of  the  non-zero  blocks 
are  associated  with  the  presence  of  all  (n!)  permutations  of  electron  assignments  to  each 
of  the  corresponding  distinct  spin-orbital  products  in  the  representation  of  Eq.  (3),  which 
circumstance  further  insures  that  any  finite  atomic  representation  in  the  form  of  Eq.  (3) 
provides  an  invariant  subspace  such  that  the  operation  of  the  antisymmetrizer  is  closed  in 
this  subspace 

PA$(r)  =  $(r)  •  PA-  (9) 

As  a  consequence,  the  eigenfunctions  of  the  matrix  representative  of  the  antisymmetrizer  in 
any  such  finite  space  will  transform  correctly  as  either  totally  antisymmetric  or  non-totally 
antisymmetric  basis  states  under  all  electron  permutations  [3]. 

The  aforementioned  eigenvalues  and  eigenfunctions  of  the  antisymmetrizer  are  obtained 
from  the  unitary  transformation  matrix  Up  that  diagonalizes  the  matrix  of  Eq.  (7), 

up  Pa  Up  =  ( (n!^tp  0°J,  (10) 

where  the  the  physical  block  (n!)  Ipp  contains  the  non-zero  eigenvalues  and  the  unphysical 
block  0UU  vanishes  identically.  The  individual  non-zero  n!-by-n!  blocks  indicated  above  are 
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each  diagonalized  separately  and  contribute  a  single  non-vanishing  eigenvalue  (n!)  each  to 
the  physical  block  of  Eq.  (10),  with  all  other  eigenvalues  zero.  The  eigenfunctions  of  the 
antisymmetrizer  take  the  correspondingly  blocked  form 

d>P(r)  =  *(r)  •  Up  =  {{*P(r)}p,  {*P(r)}u},  (11) 

where  the  physical  states  satisfy 

P.4{4>p(r)}p  =  (n!){4.P(r)}p,  (12) 

and  the  non-Pauli  states  provide  the  null  result  PJ4{$p(r)}p  =  0.  It  is  found  that  the 
physical  eigenstates  of  Eqs.  (11)  and  (12)  are  Slater  determinants  made  up  of  the  config¬ 
urations  of  spin-orbitals  comprising  the  representation  of  Eq.  (3),  whereas  the  non-Pauli 
eigenstates  are  totally  symmetric  or  are  linear  combinations  of  other  (degenerate)  repre¬ 
sentations  of  the  nth- degree  symmetric  group  [16]. 

Equations  (7)  to  (12)  indicate  that  the  spectrum  of  Pa  acting  in  the  domain  of  the  spectral- 
product  basis  of  Eq.  (3)  is  that  of  a  compact  operator,  with  zero  providing  the  lower 
limiting  point  of  the  spectrum  and  the  associated  states  {<f>p(r  :  R)}u  corresponding 
to  a  point  of  accumulation,  whereas  the  totally  antisymmetric  states  { <h p ( r  :  R)}p  are 
associated  with  the  upper  limiting  point  n!  and  correspond  to  its  point  spectrum  [17]. 

The  unitary  matrix  of  Eq.  (10)  provides  a  transformed  Hamiltonian  matrix  in  the  form 

HP,UPH.UP=(<H^  {Hp°}J,  (13) 

where  {Hp}pp  and  (Hp}uu  are  the  physical  (pp)  and  unphysical  (pm)  blocks  which  provide 
the  totally  antisymmetric  and  the  non-Pauli  solutions,  respectively,  the  off-diagonal  blocks 
vanishing  identically  in  accordance  with  Eq.  (9).  The  physical  Hamiltonian  matrix  of 
Eq.  (13)  is  identical  with  that  obtained  from  the  Slater  determinants  constructed  by 
selecting  only  linearly-independent  ordered  spin-orbital  configurations  from  the  basis  of  all 
such  products  of  Eq.  (3).  Accordingly,  the  transformation  matrix  Up  of  Eq.  (10)  obtained 
from  the  matrix  representative  of  the  antisymmetrizer  of  Eq.  (7)  correctly  incorporates  the 
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effects  of  electron  permutation  symmetry,  in  the  absence  of  prior  basis-set  antisymmetry, 
in  the  Hamiltonian  matrix  of  Eq.  (4)  through  the  transformation  of  Eq.  (13). 


2.2  Spectral-Product  Formalism  for  Polyatomic  Molecules 

The  foregoing  introductory  development  applies  with  some  important  differences  and  mod¬ 
ifications  to  the  adiabatic  (Born-Oppenheimer)  electronic  structures  of  molecules.  Specifi¬ 
cally,  all  quantities  given  in  Section  2.1  now  depend  explicitly  upon  the  spatial  arrangement 
of  atoms  in  the  molecule,  specified  collectively  by  R  =  (Ri,  R2, ... Rn ),  with  the  coordi¬ 
nates  Ra  giving  the  individual  atomic  positions.  In  order  to  follow  closely  the  atomic 
development  of  Section  2.1,  and  to  avoid  unnecessary  notational  complexity,  the  molecular 
development  is  presented  here  for  a  collection  of  n  one-electron  (hydrogen)  atoms  [9],  with 
generalization  to  many-electron  atoms  requiring  only  minor  additional  elaboration  [18]. 

The  spectral-product  representation  of  Eq.  (3)  for  polyatomic  II„  molecules  becomes 

$(r  :  R)  —  j <//■*-) (1)  0  $^(2)  0  •  •  ■  <//n)(n)j  ,  (14) 

where  the  indvidual  spectral  states  are  the  one-electron  spin-orbital  row  vectors  of 

Eq.  (3),  each  now  centered  at  the  different  origins  Ra,  with  the  coordinates  i  measured 
relative  to  these  centers.  As  in  the  atomic  case  [16],  the  basis  of  Eq.  (14)  spans  the 
totally  antisymmetric  representation  of  the  aggregate  n-electron  permutation  group  once 
and  only  once  in  the  limit  of  closure,  as  well  as  all  other  (non-Pauli)  representations  of  the 
nth -degree  symmetric  group  in  this  limit  [9]. 


The  molecular  Hamiltonian  operator  is  [cf.,  Eq.  (2)] 

n  n 

H(r  :R)  =  ^  {Hf‘-\i)+  Y  V^\v,j  :  Raf,)} 

a=l  /3=a+l 


n2  2 

— v2 

2m 


}+  E  { 


P=a+1 


Ra(3  Tj 


a  T  ij 


whereas  its  matrix  representative  in  the  basis  of  Eq.  (14)  [cf.,  Eq.  (4)], 


H(il)  =  <*(r:  R)\H(r:R)\&(r:R))  =  J2{H{a)+  E  V^ii^)}, 

a=l  /3=a+l 


(15) 


(16) 
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is  seen  to  include  atomic  terms  H*-®-*  =  (<f>(r  :  R)\H^a\i)\^(r  :  R ))  and  atomic-pair 
interaction  terms  (Rap)  =  (<fr(r  :  R)\V(-a,l3'> (i;  j  :  Rap)\&(r  :  R ))  which  provide 

Hermitian  matrix  representatives  of  the  corresponding  atomic  and  interaction  operators 
in  the  molecular  Hamiltonian  operator  of  Eq.  (15).  The  atomic  energy  matrix 
is  independent  of  atomic  position  and  includes  only  diagonal  one-electron  energies  when 
the  orthonormal  atomic  spin-orbitals  (i)  are  chosen  to  be  atomic  pseudo-eigenstates, 
whereas  the  interaction  matrix  {Rap)  is  generally  non-diagonal  and  depends  ex¬ 

plicitly  upon  the  vector  separation  Rap  of  the  two  indicated  atoms,  but  not  upon  the 
individual  laboratory-frame  positions  of  the  two  atoms,  nor  upon  the  position  vectors  of 
the  other  (n  —  2)  atoms  in  the  molecule.  The  particularly  simple  form  of  Eq.  (16)  is  largely 
a  consequence  of  the  ortho  normality  of  the  spectral-product  basis  employed  in  the  absence 
of  prior  enforcement  of  overall  electron  antisymmetry,  and  the  atomic  pairwise-additive 
nature  of  the  interaction  terms  in  the  Hamiltonian  operator  of  Eq.  (15). 

The  Hermitian  matrix  representative  of  the  antisymmetrizer  P\  constructed  in  the  basis 
of  Eq.  (14), 

P a(R)  =  (*(r  :  R)\PA\*(r  :  R))  =  (n!)-1(PA$(r  :  R)\PA®(r  :  R)),  (17) 

depends  explictly  on  the  atomic  arrangement  R.  and  is  seen  to  be  proportional  to  the 
metric  matrix  of  the  explicitly  antisymmetrized  product  basis  [13].  In  contrast  to  the 
atomic  matrix  of  Eq.  (7),  the  matrix  of  Eq.  (17)  is  a  complicated  function  of  it,  requiring 
detailed  computational  evaluation.  Moreover,  the  operation  of  the  antisymmetrizer  is 
not  closed  in  any  finite  representation  of  the  form  of  Eq.  (14),  and  requires  a  complete 
spectral  representation  to  achive  closure  P,\ 3>(r)  — ■>  <1 >(r)  •  Pa  as  a  limiting  process.  This 
circumstance  is  a  consequence  of  the  fact  that,  in  contrast  to  the  situation  in  the  atomic  case 
of  Eq.  (3),  the  representation  of  Eq.  (14)  does  not  include  explicitly  all  n\  permutations 
of  electron  assignments  to  the  spin-orbital  product  states  appearing  there,  with  exchange 
terms  included  only  implictly  in  the  limit  of  a  complete  spin-orbital  representation. 

Correspondingly,  the  unitary  transformation  matrix  Up  (it)  required  to  partition  the 
Hamiltonian  matrix  of  Eq.  (16)  into  physical  and  non-Pauli  blocks  is  obtained  in  a  limit 
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from  the  diagonalization  [cf.,  Eq.  (10)] 

UP(fl)t.pA(fl).Up(fl)=({P^>-  {Pd(°R)u)  ^  ({nf*  S).  (18) 

Here,  the  eigenvalues  of  the  matrix  Pa(-R)  are  partitioned  into  an  upper  diagonal  block 
{Pd{R)}pp  containing  the  largest  eigenvalues,  which  tend  to  the  upper  limiting  point 
(n\)  of  the  spectrum  in  the  spectral  closure  limit,  and  a  lower  diagonal  block  {P^(i2)}wu 
of  eigenvalues  which  tend  to  zero  in  this  limit  [10-12],  Similarly,  the  eigenstates  of  the 
antisymmetrizer  in  this  case  are  obtained  as  a  limit  in  the  form  [cf.,  Eq.  (12)] 

$P(r  :  R)  =  $(r  :  R)  ■  U P(R)  -  {{$P(r  :  R)}p,  {*P(r  :  i?)}w},  (19) 

where  { CIJP  (r  :  R)  }p  contains  the  totally  antisymmetric  states  corresponding  to  the  non¬ 
zero  eigenvalues  of  the  antisymmetrizer  and  {<hP(r  :  R) contains  the  non-Pauli  states 
corresponding  to  the  zero  eigenvalues  of  the  antisymmetrizer. 

Finally,  the  matrix  UP(_R)  of  Eq.  (18)  is  employed  in  constructing  the  transformed  Hamil¬ 
tonian  matrix 

n  n 

HP(^)  =  UP(i2)t-H(i?)-UP(i?)  =  ^{H^)(i?)+  ^  V^’^(fl)} 

a=l  /3=a+l 

__  ( {Hp(i?)}pp  {Hp(H)}p„\  /{Hp(H)}pp  0  \ 

{Hp(iC)}„u  J  (  0  {HP  (fl)}„„p  ’ 

where  Hp  '  (R)  and  Vp  ’^(R)  refer  to  transformations  of  the  individual  atomic  and  in¬ 
teraction  matrices  in  Eq.  (16),  {HP(i?)}pp  and  {HP(i?)}uu  are  the  physical  (pp)  and 
unphysical  (mm)  blocks  of  the  transformed  Hamiltonian  matrix  which  provide  the  antisym¬ 
metric  and  the  non-Pauli  solutions,  respectively,  and  the  off-diagonal  blocks  vanish  in  the 
limit  (— >)  of  spectral  closure  [9-12],  The  transformation  of  Eq.  (20)  is  seen  to  incorpo¬ 
rate  the  non-local  effects  of  overall  electron  antisymmetry  on  the  atomic  H^b  — >  Hp  "*  (R) 
and  interaction  \^a,P\Rap)  — >  V^*’^(.R)  matrices,  which  consequently  now  depend  on 
the  positions  of  all  atoms  in  the  molecule.  Nevertheless,  the  Hermitian  atomic  and  pair- 
interaction  matrices  of  Eq.  (20)  individually  have  appropriate  dissociation  limits,  with  the 
entire  Hamiltonian  matrix  approaching  the  sum  of  the  individual  atomic  terms  of  Eq.  (16) 
in  the  complete  dissociation  limit  R  — >  oo. 
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2.3  Atomic-Pair  Implementation 

The  particularly  simple  forms  of  the  atomic  and  interaction-energy  terms  appearing  in  the 
Hamiltonian  matrix  of  Eq.  (16)  suggest  the  possibility  of  their  systematic  evaluation  in  a 
series  of  individual  atomic  and  atomic-pair  calculations  which  can  be  performed  once  and 
for  all  independently  of  any  particular  atomic  aggregate  arrangement.  Since  atomic  and 
diatomic  eigenspectra  are  generally  available  from  experimental  and  previous  theoretical 
studies,  it  is  convenient  to  cast  such  an  approach  in  the  context  of  calculations  of  atomic  and 
diatomic  eigenstates  and  energies,  which  can  be  critically  evaluated  by  suitable  comparisons 
prior  to  their  incorporation  in  the  spectral- product  representation  of  Eqs.  (14)  to  (20). 

To  facilitate  the  introduction  of  fragment  information,  the  Hamiltonian  operator  of  Eq. 
(15)  and  correponding  matrix  representative  of  Eq.  (16)  can  be  written  in  the  alternative 
but  equivalent  forms  [6] 

n  n 

H(r  :R)  =  J2  ^  :  Rap)  -  H^\i)  -  (21) 

a=l  /3=a+l 

and 

n  n 

H(R)  =  {  H(q)  +  {K{a,l3){Rap)  -H(a)  -H(/3)}|,  (22) 

a=l  /3=a+ 1 

where  :  Rap)  =  H^a\i)  +  H^\j)  +  :  Rap)  is  the  diatomic  Hamil¬ 
tonian  operator  for  atoms  a  and  /?,  and  H (a’f3') (Rap)  —  (Rap)  is 

its  corresponding  matrix  representative  in  the  basis  of  Eq.  (14).  These  expressions  are 
sufficient  to  relate  calculations  on  individual  atoms  and  diatomic  molecules  in  the  repre¬ 
sentation  of  Eq.  (14)  to  the  polyatomic  Hamiltonian  matrix  of  Eqs.  (16)  and  (20)  also 
constructed  in  this  basis. 


The  physical  block  of  the  aggregate  Hamiltonian  matrix  of  Eq.  (20)  constructed  in  this 
manner  takes  the  form  [18] 

n  n 

{HP(B)}OT  =  ^{{H‘?)(iJ)}OT+  ^  {H^'3)(K)-H^,(Ji)-H‘f)(fl)}i>p}.  (23) 

a=l  /3=a+l 

where  the  atomic  matrices  {Hp^(ii)}pp  and  {H jf}(ii)}  are  defined  in  connection  with 
Eq.  (20)  above,  and  the  physical  block  of  the  diatomic  Hamiltonian  matrix  is 

{H<,“'»(iJ)}„,  =  {U^Jl)}^  ■  {H<?-«(^)}pP  •  (24) 
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Evaluation  of  Eq.  (24)  requires  the  transformation  matrix 


{U<?'W («)}„,  =  {(4“J3)V,<  :  «)l*p(r  :  R))}„,  (25) 

where  \rd  :  R) } p  and  {<l>p(r  :  R)}p  are  the  totally  antisymmetric  orthonormal 

states  of  Eq.  (19)  for  the  diatomic  molecule  (a,  (3)  and  for  the  polyatomic  molecule  of 
interest,  respectively,  whereas  the  diatomic  Hamiltonian  matrix  in  Eq.  (24)  is  given  by 

{H =  {U {£m(Raf,)}m  ■  •  {U (26) 

where  {E  ^\RaP)}pp  is  the  diatomic  eigenvalue  matrix  constructed  in  the  representa¬ 
tion  of  Eq.  (14)  and  {U^ ’^(i?a/3)}pp  is  the  unitary  matrix  which  provides  the  diatomic 
eigenfunctions  in  the  form  {^r(a>/3)(rd  :  R)}p  —  {<&{p,l3\rd  :  R)}p  ■  {U (Ra/3)}pp-  The 
diatomic  states  are  functions  of  electron  coordinates  ( rd )  in  the  diatomic  frame,  as  in¬ 
dicated,  and  all  matrices  and  vectors  in  Eqs.  (23)  to  (26)  have  the  dimensions  of  the 
representation  of  Eq.  (14),  although,  with  the  exception  of  the  polyatomic  Pa(-R)  matrix 
of  Eq.  (17),  their  evaluations  require  only  explicitly  atomic  and  diatomic  calculations  [18]. 

The  diatomic  Hamiltonian  matrix  of  Eq.  (24)  is  seen  to  include  modifications  of  the  cor¬ 
responding  bare  diatomic  energy  matrix  appearing  in  Eq.  (26)  due  to  the  inter-atomic 
antisymmetry  required  upon  of  immersion  of  the  interacting  atomic  pair  in  the  molecular 
aggregate,  effects  which  are  incorporated  in  the  development  through  the  indicated  trans¬ 
formation  matrix.  Accordingly,  Eqs.  (23)  and  (24)  can  provide  an  exact  representation 
of  the  entire  polyatomic  Hamiltonian  matrix,  in  the  limit  of  closure,  in  the  form  of  sums 
of  Hermit ian  atomic  and  atomic-pair-interaction  matrices  which  can  be  determined  from 
atomic  and  diatomic  calculations  peformed  once  and  for  all  independently  of  any  particular 
atomic  arrangement. 

3.0  Illustrative  Computational  Applications 

Illustrative  computational  applications  are  made  of  the  atomic-pair  implementation  de¬ 
scribed  in  the  preceeding  section  in  the  case  of  the  H3  molecule  employing  a  minimal 
Is3  basis  [19]  and  a  12-term  Gaussian  representation  of  the  Is  atomic  hydrogen  orbital. 
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This  basis  provides  an  orbital  energy  (-0.4999  au)  in  good  agreement  with  the  correct 
value,  as  well  as  explicitly  antisymmetric  singlet  and  triplet  molecular  hydrogen  eigenstates 
{\l/ {a,P)(rd  :  R)}p  and  a  ground-state  equilibrium  interatomic  separation  ( Re  =  1.642  a o) 
and  total  energy  (-1.1157  au)  in  good  agreement  with  the  accepted  Heitler-London  values. 
Conventional  Gaussian-based  computational  methods  are  employed  in  evaluations  of  all 
quantities  required  in  this  application  of  the  atomic-pair  development  for  H3  [2], 

In  Fig.  1  are  shown  contour  plots  of  the  two  doublet  eigenvalues  of  the  Pa  (R)  matrix 
of  Eq.  (17)  obtained  in  the  Is3  H3  representation  for  non-symmetric  collinear  atomic 
arrangements.  The  eigenvalue  surface  of  Fig.  1(a)  has  a  maximum  near  the  saddle  point 
( Rab  —  Rac  —  1-80  ao),  decreases  rapidly  to  zero  in  the  limit  Rai,  =  Rac  — »  0,  approaches 
1  in  the  three-body  break-up  limit  Rab  —  Rac  — >  00,  and  goes  to  2  in  the  diatomic 
limits  Rab  — >  0,  Rac  — >  00  and  Rab  — >  00 ,  Rac  — 1 ►  0.  The  vanishing  of  the  eigenvalue  of 
Fig.  1(a)  as  Rab  and  Rac  — »  0,  corresponding  to  formation  of  an  invalid  Is3  atomic  Li 
configuration,  indicates  that  the  minimal  Is3  basis  does  not  support  a  physically  significant 
H3  state  in  this  limit.  The  smaller  doublet  eigenvalue  surface  shown  in  Fig.  1(b)  is 
evidently  less  structured  than  that  of  Fig.  1  (a)  and  vanishes  in  both  the  diatomic  and 
triatomic  united-atom  limits.  The  Pa  (R)  -matrix  eigenvalue  surfaces  evidently  identify 
atomic  configurations  in  which  improper  states  arise,  provide  information  complementary 
to  the  associated  energy  surfaces,  and  can  anticipate  the  presence  of  crossings  in  the 
corresponding  potential  energy  surfaces,  as  demonstrated  in  further  detail  below. 

Figure  2  provides  constant-energy  contour  plots  of  the  two  lowest-lying  doublet  energy 
surfaces  in  H3  obtained  from  Eqs.  (23)  to  (26)  in  the  Is3  representation  for  the  non- 
symmetric  collinear  atomic  arrangements  of  Fig.  1.  These  well-known  energy  surfaces 
are  presented  here  to  illustrate  the  extent  of  agreement  obtained  with  corresponding  con¬ 
ventional  valence-bond  calculations  in  a  Is3  basis  [20],  and  to  demonstrate  their  general 
accordance  with  accurate  quantum-chemistry  calculations  [21-24],  Specifically,  the  calcu¬ 
lated  saddle-point  energy  (-1.574  au)  of  Fig.  2(a)  compares  well  with  that  obtained  from 
valence-bond  calculations  (-1.583  au)  in  the  same  basis,  and  the  H3  geometry  at  the  saddle 
point  ( Rab  =  Rac  —  2.18  ao)  is  in  general  accord  with  the  accepted  position  (Rab  —  Rac 
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=  1.757  ao)  [21-24],  The  current  calculation  places  the  saddle  point  approximately  0.0735 
au  (2.00  eV)  below  the  three-body  (H+H+H)  breakup  energy  of  -1.50  au.  and  thus  pro¬ 
vides  approximately  50%  of  the  accepted  accurate  saddle-point  binding  energy  of  0.1570 
au  (4.273  eV)  relative  to  this  limit  [21-24],  The  calculated  saddle-point  energy  relative  to 
the  two-body  break-up  energy  (H-f  II2  j  provides  a  1.14  eV  barrier  to  the  exchange  reaction 
which  is  in  accord  with  the  valence-bond  calculations  (0.876  eV)  but  is  about  three  times 
larger  than  the  accepted  value  of  0.4166  eV  [21-24],  emphasizing  the  well-known  sensitivity 
of  this  energy  difference  to  the  computational  approximation  employed. 

The  structure  of  the  ground-state  energy  surface  of  Fig.  2(a)  is  evidently  complementary 
to  that  of  the  associated  PA(-R)-matrix  eigenvalue  surface  of  Fig.  1(a),  with  larger  values 
of  the  Pa(-R)  -matrix  eigenvalues  corresponding  to  generally  lower  energy  values.  Sim¬ 
ilarly,  the  energy  surface  for  the  excited  doublet  state  shown  in  Fig.  2(b)  is  seen  to  be 
monotonically  increasing  from  the  three-atom  break-up  limit  as  the  diatomic  and  tri atomic 
united-atom  limits  are  approached,  in  accord  with  the  monotonically  decreasing  form  of 
the  associated  Pa(-R)  -matrix  eigenvalue  surface  of  Fig.  1(b). 

In  Fig.  3  are  shown  comparisons  of  the  first  two  doublet  energies  in  H3  with  correspond¬ 
ing  valence-bond  calculations  for  symmetric  collinear  (Rab  —  Rac )  atomic  arrangements. 
Evidently,  the  present  results  for  the  ground-state  energy  surface  are  in  good  agreement 
with  the  Is3  valence-bond  calculations  which  include  only  covalent  structures  (o),  and  are 
also  in  accord  with  valence-bond  results  which  include  all  possible  ionic  configurations  (•) 
in  the  Is3  basis  [20],  The  first  excited  doublet  state  in  H3  is  seen  to  be  repulsive,  in  good 
agreement  with  the  valence-bond  calculations  and  in  qualitative  accord  with  previously  re¬ 
ported  accurate  calculations,  although  this  excited  doublet  is  known  to  undergo  an  avoided 
crossing  with  a  higher-lying  state  which  is  not  included  in  the  present  development  [23]. 
Although  the  results  of  Fig.  3  for  the  ground-state  energy  surface  appear  to  bound  the 
valence-bond  results  from  above,  this  is  only  a  happenstance  in  view  of  the  small  basis 
set  employed.  Indeed,  similar  apparent  boundedness  is  not  found  for  the  excited  doublet 
state  in  Fig.  3,  although  both  results  are  in  satisfactory  agreement  with  the  valence-bond 
calculations,  the  ionic  configurations  making  generally  modest  contributions  to  both  states. 
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Figure  4  depicts  contour  plots  of  the  two  largest  doublet-state  Pa  (-R)- matrix  eigenvalues 
in  H3  for  C-2v  arrangements  in  which  the  Ra},  —  Rac  separations  are  co-varied  and  the 
apex  angle  6  is  varied  from  30  to  180  degrees.  A  seam  of  intersection  of  these  two  surfaces 
indicated  by  the  cusps  in  the  contours  is  evident  in  the  figures  at  6  =  60  degrees,  which 
corresponds  to  the  familiar  degeneracy  associated  with  D^h  symmetry  [25].  The  dashed 
and  solid  portions  of  the  contours  in  the  two  figures  depict  the  continuous  surfaces  which 
pass  smoothly  through  one  another  at  the  6  —  60  degrees  line  of  intersection.  In  panel  (a) 
the  higher-value  eigensurface  (solid  lines)  is  to  the  right  of  the  60  degree  line  and  the  lower- 
value  eigensurface  (dashed  lines)  to  the  left  of  this  line  of  intersection.  Correspondingly, 
panel  (b)  shows  the  reversed  arrangement,  with  the  lower  eigenvalue  surface  (dashed  lines) 
to  the  right  and  the  higher  (solid  lines)  to  the  left  of  the  60  degree  line.  The  associated 
energy  surfaces  obtained  from  Eqs.  (23)  to  (26)  for  the  ground-  and  first-excited  doublet 
states  of  H3  for  C^v  atomic  arrangements,  shown  in  Fig.  5,  similarly  depict  the  expected 
D?,h  degeneracy  at  6  =  60  degrees  and  the  associated  surface  crossing  along  the  cusps  in 
the  contours,  complementary  to  the  corresponding  seam  of  intersection  anticipated  by  the 
Pa(A)  -matrix  eigenvalues  of  Fig.  4. 

In  Fig.  6  are  shown  comparisons  of  the  two  lowest-lying  doublet  energies  in  H3  with 
coresponding  valence-bond  calculations  for  the  C^v  atomic  arrangements  of  Fig.  5  in 
which  Raf,  and  Rac  are  held  fixed  at  1.64  a 0  and  6  is  varied  from  20  to  180  degrees. 
Evidently,  the  present  results  are  in  good  agreement  with  the  corresponding  Is3  valence- 
bond  calculations  which  include  only  covalent  structures  (o).  By  contrast,  although  the 
present  results  are  also  in  accord  with  valence-bond  results  which  include  all  possible  ionic 
configurations  (•)  in  the  Is3  basis  [20],  the  undulation  between  120  and  180  degrees  evident 
upon  inclusion  of  such  terms  in  the  valence-bond  calculations  is  evidently  not  present  in  the 
Is3  implementation  of  the  atomic-pair  formalism.  The  crossing  of  the  two  surfaces  when  0 
=  60  degrees  is  evident  in  the  figure,  and  in  accord  with  the  valence-bond  calculations.  As 
in  Fig.  3,  the  apparent  upper  bound  provided  by  the  present  results  on  the  valence-bond 
calculations  in  the  case  of  the  ground  state  in  Fig.  6  is  a  happenstance  in  view  of  the  small 
basis  representation  employed,  with  a  similar  apparent  upper  boundedness  not  present  in 
the  case  of  the  excited  doublet  state. 
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4.0  Concluding  Remarks 


The  theoretical  development  reported  in  Section  2  provides  an  alternative  perspective  on 
electronic  structure  calculations  for  atoms  and  molecules  and  other  forms  of  matter.  Anti¬ 
symmetry  restrictions  are  enforced  subsequent  to  construction  of  the  many-electron  Hamil¬ 
tonian  matrix  for  an  atom  or  molecule  in  an  orthogonal  spectral-product  basis.  Transfor¬ 
mation  to  a  permutation-symmetry  representation  obtained  from  the  eigenstates  of  the 
aggregate  electron  antisymmetrizer  enforces  the  requirements  of  the  Pauli  principle  in  this 
approach.  In  this  way,  results  are  obtained  in  applications  to  many-electron  atoms  which 
are  identical  with  the  use  of  antisymmetrized  configurational  state  functions  in  variational 
calculations,  providing  some  degree  of  confidence  in  the  soundness  of  the  method  more 
generally.  In  applications  to  polyatomic  molecules,  the  development  accommodates  the 
incorporation  of  fragment  information  in  the  form  of  Hermitian  matrix  representatives 
of  atomic  and  diatomic  operators  which  include  the  non-local  effects  of  overall  electron 
antisymmetry  on  the  individual  atomic  and  pairwise- atomic  interaction  terms. 

Illustrative  computational  applications  of  the  development  to  the  H3  molecule  reported 
in  Section  3  employ  largely  conventional  methods  of  calculations  for  atoms  and  diatomic 
molecules,  and  can  benefit  from  the  available  highly-developed  methods  devised  for  this 
purpose.  The  eigensurfaces  of  the  matrix  representative  of  the  electron  antisymmetrizer 
are  seen  to  identify  the  presence  of  non-Pauli  states  in  the  spectral-product  basis  employed, 
to  complement  the  structures  of  the  more  familiar  energy  surfaces,  and  to  anticipate  the 
presence  of  seams  of  surface  intersections  generally  associated  with  high-symmetry  molec¬ 
ular  geometries.  The  minimal-basis-set  results  reported  provide  energy  surfaces  in  good 
agreement  with  corresponding  valence-bond  calculations,  and  are  in  general  accord  with 
the  results  of  accurate  quantum-chemistry  calculations. 

The  formalism  and  applications  reported  here  avoid  repeated  calculations  of  the  one-  and 
two-electron  integrals  generally  required  in  construction  of  polyatomic  Hamiltonian  matri¬ 
ces  employing  the  antisymmetric  basis  states  commonly  adopted  in  calculations  of  potential 
energy  surfaces.  Accordingly,  the  approach  can  possibly  provide  an  alternative  ab  initio 
method  suitable  for  computational  applications  to  polyatomic  molecules  more  generally. 
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Figure  1-Constant-value  contours  for  the  two  largest  doublet  eigenvalues  of  the  Pa(-R) 
matrix  of  Eq.  (17)  for  collinear  atomic  configurations  of  H3  constructed  in  a  Is3  rep¬ 
resentation,  employing  an  increment /decrement  of  0.05  between  adjacent  contours.  The 
coordinate  axes  give  the  distances  Rab  and  Rac  between  the  central  atom  located  at  the 
origin  and  the  outer  atoms,  measured  in  atomic  units  (ao).  In  panel  (a)  the  contours 
increase  uniformly  in  value  from  that  labeled  1.05,  except  as  R(li,  =  Rac  — >  0  inside  the 
saddle  point,  in  which  region  the  values  decrease  monotonically  to  zero,  whereas  in  panel 
(b)  the  contour  values  decrease  uniformly  from  that  labeled  0.95. 
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Figure  2-Constant-energy  contours  for  the  ground  and  first  excited  doublet-state  poten¬ 
tial  energy  surfaces  of  H3  for  non-symmetric  collinear  atomic  arrangements,  obtained  from 
Eqs.  (23)  to  (26)  in  a  Is3  representation,  employing  an  energy  decrement /increment  of 
0.5  eV  between  contours.  The  coordinate  axes  are  as  in  Fig.  1,  and  the  zero  of  energy 
employed  (-1.50  cm)  is  that  of  the  three-atom  dissociation  limt  (H+H+H).  Panel  (a)  de¬ 
picts  the  familiar  saddle  and  diatomic  structures  of  the  ground-state  energy  surface  in  H3 
[25],  whereas  the  excited-state  energy  surface  of  panel  (b)  is  montonically  increasing  with 
decreasing  separations  Rab  and/or  Rac. 
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Figure  3-Ground  and  first-excited  doublet  state  energies  in  H3  for  symmetric  collinear 
atomic  arrangements.  The  solid  curves  depict  the  values  obtained  from  Eqs.  (23)  to  (26) 
in  a  Is3  representation,  whereas  the  data  points  are  obtained  from  conventional  valence- 
bond  calculations  in  a  Is3  basis  set  including  (•)  and  excluding  (o)  ionic  configurations 
[20].  The  zero  of  energy  employed  (-1.50  au)  is  that  of  the  three-atom  dissociation  limit 
(H+H+H),  as  in  Fig.  2. 
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Figure  4-Constant-value  contours  for  the  two  largest  doublet-state  eigenvalues  of  the 
Pa(-R)  matrix  of  Eq.  (17)  in  a  Is3  representation  for  C^v  atomic  arrangements  (Rab  = 
Rac),  where  0  is  the  apex  angle  of  the  tri atomic  H3  configuration,  employing  an  incre¬ 
ment/decrement  of  0.1  between  adjacent  contours.  The  cusps  in  the  contours  of  both 
panels  indicate  the  presence  of  a  seam  of  intersection  at  6  =  60  degrees  in  two  continu¬ 
ous  surfaces  identified  by  the  solid  and  dashed  lines,  with  the  contour  values  decreasing 
uniformly  from  right  to  left  in  panel  (a)  and  decreasing  uniformly  from  top  to  bottom  in 
panel  (b),  as  is  discussed  further  in  the  text. 
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Figure  5-Constant-energy  contours  for  the  two  lowest-lying  doublet  states  of  H3  for 
C2V  atomic  arrangements  (Rab  =  Rac ),  obtained  from  Eqs.  (23)  to  (26)  in  a  Is3  repre¬ 
sentation,  where  6  is  the  apex  angle  of  the  tri atomic  configuration,  employing  an  energy 
increment /decrement  of  0.5  eV  between  adjacent  contours.  The  zero  of  energy  (-1.50  au ) 
is  that  of  the  three-atom  dissociation  limit  (H+H+H).  The  cusps  in  the  contours  of  both 
panels  indicate  the  familiar  D^h  symmetry  seam  of  intersection  at  9  =  60  degrees  in  the  two 
continuous  energy  surfaces  identified  by  the  solid  and  dashed  lines  [25],  with  the  contour 
energy  values  increasing  uniformly  from  top  to  bottom  in  both  panels.  Panel  (a)  depicts 
the  lowest-lying  surface  only  for  values  of  9  >  60  degrees,  whereas  panel  (b)  depicts  the 
lowest-lying  surface  only  for  6  <  60  degrees,  as  is  discussed  further  in  the  text. 
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Figure  6-Ground  and  first-excited  doublet-state  energies  in  H3  for  C^v  atomic  arrange¬ 
ments  (Rab  —  Rac  —  1-64  ao),  where  6  is  the  apex  angle  of  the  triatomic  configuration. 
The  solid  curves  depict  the  values  obtained  from  Eqs.  (23)  to  (26)  in  a  Is3  representation, 
whereas  the  data  points  are  obtained  from  conventional  valence-bond  calculations  in  a  Is3 
basis  set  including  (•)  and  excluding  (o)  ionic  configurations  [20].  The  zero  of  energy  (-1.50 
au)  is  that  of  the  three-atom  dissociation  limit  (H+H+H),  as  in  Fig.  5.  The  crossing  point 
of  the  two  curves  corresponds  to  D3h  symmetry,  as  discussed  further  in  the  text. 
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